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' The Bianchi IX model has been used often to investigate the structure close to singularities 

^ . of general relativity. Its classical chaos is expected to have, via the BKL scenario, implications 

\ even for the approach to general inhomogeneous singularities. Thus, it is a popular model to test 
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consequences of modifications to general relativity suggested by quantum theories of gravity. This 
^ ■ paper presents a detailed proof that modifications coming from loop quantum gravity lead to a 
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. non-chaotic effective behavior. The way this is realized, independently of quantization ambiguities, 
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\ suggests a new look at initial and final singularities. 
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I. INTRODUCTION 

It is well known that according to the singularity theorems Q], the space-time describ- 
ing the backward evolution of an expanding universe is necessarily singular in the sense of 
geodesic incompleteness. The generality of singularity theorems also prevents them from 
giving any information about the nature of such a singularity for example in terms of cur- 
vature invariants. By contrast, the BKL approach to the issue of singularities asks [2]: In a 
neighborhood of a presumed singularity, is there a general solution of the Einstein equation 
such that at least some curvature invariants diverge? In this formulation of the question, one 
also obtains information about the possible nature of the presumed singularity described in 
terms of the approach to the singularity. The conclusion of the long and detailed analysis is 
summarized in the BKL scenario: Generically, as the singularity is approached, the spatial 
geometry can be viewed as a collection of small patches each of which evolves essentially 
independently of the others according to the Bianchi IX evolution (for recent numerical 
evidence see jj^l). The Bianchi IX evolution towards its singularity is described by an infi- 
nite succession of Kasner evolutions (two directions contracting while third one expanding) 
punctuated by permutations of the expanding/contracting directions as well as possible ro- 
tations of the three directions themselves. The qualitative analysis of these permutations 
and rotations ('oscillations') is done very conveniently in terms of a billiard ball bouncing 
off moving walls and has been analyzed for a possible chaotic behavior ^. 

However, it has been well recognized that the conclusion of ad infinitum oscillations of 
the BKL singularity, a consequence of unbounded growth of the spatial curvature, cannot 
be trusted very close to the singularity where the classical Einstein equations themselves are 
expected to break down. Presumably the equations will be superceeded by some quantum 
extensions. In the absence of any specific and detailed enough quantum theory the questions 
of the fate of the classical singularity vis a vis the behavior of space-time near such a 
region, could not have been addressed. The focus therefore has been to include qualitatively 
expected quantum modification. For example, in the Kaluza-Klein picture and more recently 
the stringy picture, the minimal expected modification is the matter content - notably the 
dilaton and p-form fields. Higher derivative corrections in the effective action are also 
expected. Another qualitative implication of string theory namely the brane world scenarios 
is also a candidate to study implications of modified Einstein equations. All of these have 
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been explored with varying conclusions regarding the BKL behavior, e.g. 0,0]. However, 
they still contain the singularity and hence all these modified equations must break down 
near the singularity raising questions about the validity of the conclusion. 

In recent times this situation has changed, particularly in the cosmo logical context. 
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uce tne nature oi quantum enects. uue to tne unaeriymg quantum geometry 
Q, the internal time (volume) evolution becomes a difference equation HQ 
and one also has a quantization of the inverse scale factor with a bounded spectrum 
Thus the classically indicated singularity is removed i.e. the evolution extends to 'negative' 
(internal) time and the spatial curvatures remain finite 16, 17, 18, 3|- The equations also 



admit a continuum approximation thereby recovering the classical continuum geometry for 
larger volumes j20|, |21 1 . 

It is now meaningful to ask if and how the BKL behavior is modified. There is in fact 
a hint of what may be expected. For larger volumes, one can trust the classical picture. 
Then, in the general inhomogeneous cosmo logical context, one can approximate smaller 
patches of the spatial geometry by the homogeneous Bianchi IX models. As the volume 
is decreased, the patches have to be made smaller to sustain the approximation. // the 
BKL behavior were to continue for these individual patches, then the fragmentation into 
smaller patches must continue ad infinitum. But the underlying discrete structure cannot 
support such infinite fragmentation. Therefore, the quantum geometry which is responsible 
for singularity removal, must also ensure that Bianchi IX behavior may at the most have 
finitely many oscillations. This is in fact a self consistency requirement on the procedures 
of LQC. As reported in j^], the self consistency indeed holds and we elaborate the proof in 
this paper. 

The paper is organized as follows. 

In section II, we recall the relevant details of the classical and the loop quantized descrip- 
tion of the vacuum Bianchi IX model. We describe the kinematical set-up, the Hamiltonian 
constraint and the definition of inverse triad components including a class of quantization 
ambiguities. 

In section III, we recall the continuum approximation leading to the Wheeler-DeWitt 
differential equation. We describe how the effective continuum Hamiltonian is obtained and 
how the most important non-perturbative corrections from the underlying quantum equations 
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are incorporated. We also comment on the validity of the effective description. 

In section IV, we display the effective Hamiltonian for the Bianchi IX model and specify 
the dynamics with respect to an internal time obtained by solving the effective Hamiltonian 
constraint. 

In section V, we analyze the potential term in the Hamiltonian for internal dynamics. 
This section is divided into four subsections. In subsection A, we describe the 'walls' (pos- 
itive potential) of the potential. Unhke the classical one, this modified potential has a 
non-factorizable dependence on internal time. The walls therefore not only move but also 
change their extent and height, eventually disappearing completely. Subsection B contains 
a demonstration of the stability of the Kasncr solution when the effective potential is turned 
on. In subsection C, we relate the Kasncr solution in terms of the internal time to the more 
common description in terms of the scale factor evolution. In the last subsection, we specify 
the domain of validity of the effective description in quantitative terms. 

In section VI, we extrapolate beyond the domain of validity mainly to contrast the usual 
oscillatory behavior under the influence of the classical potential. The (small volume) asymp- 
totic analysis of this section uses qualitative arguments and proceeds in two main steps. In 
the first step, we rule out motion in the anisotropy plane being confined to a bounded re- 
gion. In the second step we show that the asymptotic trajectory must approach a Kasner 
trajectory exponentially. 

In section VII, we summarize our results. We also point out the insensitivity of the results 
to quantization ambiguities and discuss possible implications of the modified dynamics as 
to the structure of initial and final singularities. 

II. THE BIANCHI IX MODEL 

Bianchi models are a subclass of spatially homogeneous models such that the symmetry 
group contains a subgroup with simply transitive action on the spatial manifold. The simply 
transitive subgroups are classified in terms of three integers, n^, parameterizing the structure 
constants as 

Bianchi IX class of models have — 1, 1 — 1,2,3. Clearly, Cjj — implying that 
Bianchi IX is a class A model admitting a canonical formulation. 
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Given the subgroup acting simply transitively, one introduces the corresponding left- 
invariant 1-forms Lij^ — a;^dx" and their dual, density weighted vector fields Xj — Xjda, on 
the spatial manifold E. The 1-forms satisfy the Maurer-Cartan equations. 

In the metric description, the metrics of Bianchi IX models have the form, 

ds^ = dt'^ --fij{t)(jJ^CU-^, (1) 

while in the connection formulation, the geometry is described in terms of the invariant 
connections and densitized triads given by, 

with spatially constant and canonically conjugate 0}, pf. 

Note that while the symmetry group 'fixes' the invariant 1-forms only up to orthogonal 
transformations which preserve the structure constants of Bianchi IX, this freedom can be 
used only at one instant of time during an evolution. The t dependence resides in the 
symmetric 7/j in the metric formulation and in (f)],pf in the connection formulation. 

In the vacuum case, the Einstein equations Roa — 0, imply that the matrix 777 can be 
taken to be diagonal so that only three degrees of freedom remain in the metric formulation. 
In the presence of matter, the metric is in general non-diagonal. However, as argued by 
BKL, near the singularity, the non-diagonal terms do not affect either the nature of Kasner 
epochs (time intervals during which scales of spatial distances along three specific directions 
change as in the case of the Kasner solution) , or the 'law of replacement of Kasner exponents' . 
Instead, these terms lead to rotation of the Kasner axes themselves during a change over from 
one Kasner epoch to another. Thus, to study the 'oscillatory behavior' implied by succession 
of Kasner epochs with changes in the exponents, it is sufficient to restrict to the sub-class 
of diagonal Bianchi IX, and anyway for the vacuum case this restriction is automatically 
implied. One now has only three degrees of freedom in the metric formulation. 

In the connection formulation the situation is analogous. To begin with one has nine 
degrees of freedom, 0} but with a three parameter gauge freedom of SU(2) rotations of the 
index i. A further freedom, depending upon the parameters , would have been available if 
the left-invariant 1-forms were not thought of as a background structure in parameterizing 
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the invariant connections. To get the same three degrees of freedom as in the diagonal metric 
formulation, one then restricts the (f)\ to a "diagonal form", 0} := C(/)A} and correspondingly 
pI := p^A^^ where the S0(3)-matrix A includes gauge degrees of freedom (see for 
details). Such a restriction specifies the diagonal models in the connection formulation. 
This is not a symmetry reduction in the same sense as restriction to isotropic models is. 
Nevertheless, these restricted models can be analyzed by following procedures similar to 
those in the context of symmetry reductions [2^. As in reductions to isotropic models, this 
leads to significant simphfications in such diagonalized models because the gauge parameters 
in A and (almost; see below) gauge invariant parameters, c/, are neatly separated [ffij ]. 

A. Classical Framework 

The basic variables for diagonal Bianchi Class A models are specified via 

= c^j^I^jui (2) 

in terms of the 'gauge invariant' coefficients c/ with the pure gauge degrees of freedom 
contained in the S'0(3)-matrix A. The components c/ are not completely gauge invariant 
but subject to residual discrete gauge transformations which change the sign of two of the 
three components simultaneously. A diagonal densitized triad has the form 

= p^'^AlX^ (3) 

where Xj are left-invariant densitized vector fields dual to . Being an S0(3)-matrix, A 
satisfies 

A}Af = (5/ , eijkA^jA^'jA^K = ^uk 

The triad components p^ are subject to the same residual gauge transformations as the 
connection components c/ to which they are conjugate with basic Poisson bracket 

{c,,/} = 7«:5/ (4) 

where 7 is the Barbero-Immirzi parameter and k, = 8ttG the gravitational constant. We 
will use the value 7 = ~ 0.13 fixed by the black hole entropy calculations j^iHl- 
A diagonal co-triad has the form ejj = a(/)A}co'f with 

= 10203! sgn(ai) , = [aiaal sgn(a2) , = |aia2| sgn(a3) . (5) 



Note that the components as well as aj can take negative values, but only the overall sign 
sgn(p^p^p^) (i.e., the orientation) and the absolute values are gauge invariant. With this 
form of the co-triad, we obtain, in fact, a diagonal metric 
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The extrinsic curvature and the spin connection components likewise have a diagonal 
form: fc* := _ft'(/)A}co'f and := r(/)A}c<j^ with Kf = —\ai which appear in the relation 
C/ = Tj—'-yKi defining the Ashtekar connection components C/ in terms of the spin connection 
components F/ and the extrinsic curvature. The spin connection components for Bianchi 
IX are given by 0|, 



aj 
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ax 


aj 


ax 


aj 




pK 


+ 


pJ 


pJpK 


pj 







(6) 

for (/, J, K) an even permutation of (1, 2, 3). (7) 



Note that in contrast to the full theory, the spin connection is a covariant object within a 
homogeneous model since coordinate transformations have to respect the symmetry. Con- 
sequently, if non-zero, it cannot be made small by choosing appropriate local coordinates. 
The Hamiltonian constraint 16] is given by, 

H = 2k-^ { [(csFs + CsFs - r2r3)(l + 7"') - ci - T^c^c^] ai 

+ [(CiFg + CgFi - rir3)(l + 7"^) - C2 - l^'^CiC^] 02 

+ [(ciFa + - riFa) (1 + 7"') - C3 - ^-^CiC^] ] (8) 
= 2k'^ [(raFs - ri)ai + iViT^ - 12)02 + - T^)a^ 

-1(01(2203 + 020103 + 030.10-2)] • (9) 

In order to derive the classical field equations (for which we consider only positive o/ and 
p^) it is advantageous to transform to new canonical variables 

TT/ := 2Kjp^^'' = — a/ttj ^ai0203 = —(logo/)' and 

:= |logp^ such that {q^,nj} = K5j (10) 

where the prime denotes a derivative with respect to a new time coordinate r related to t by 
dt = ai0203dr (corresponding to the lapse function N = 010203). With these new variables 
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we have {vr/,p"^} = —2Kp^^^Sj and 

kNH = Ka,a2asH = 2 [pYiT^T^ - T^)+py{T,Ts - T2)+py{T2Ts - Ti) 

-i(7ri7r2 + vriVTs + 7r27r3)] (11) 



P P \ I P P \ I P P 

pi J \ p^ ) \ p^ 



{P'f - {P'f - {P'f - \{^l^2 + VTiTTs + TT^VTs) (12) 



1/4,4, 4\ 22 22 221/ , , \ f-\r>\ 

= 2 (C^l + ^^2 + %j ~ '^l'^2 ~ '^l'^3 '^2% ~ 2y^l'^2 + T^lT^'i + T^lT^V ■ (l>jj 

Now one can easily derive the equations of motion, 

(log a/)" = -TT^ -{'Ki,aia2a^H} 

which yield 



(14) 





[ai)" = 


(al 


-air 


4 

— a^ 




;«2)" = 


{al 


-air 


4 

- ^2 




I as)" = 


{al 


-air 


4 



or, using (q^)" = —2{'^'2 + 



37) 



(gO" = -4ayr,. 



For the Bianchi I model we would have the right hand sides of equation (|T3j) being zero 
so that log a/ = a/(r — tqj). This implies the Kasner behavior aj oc f^' where t = and 
the constraint feq. (fTT|) with Vj = 0) requires = 0:102 + OiOs + Q;2tt3 = |((cn + 02 + cts)^ — 
af — — a|). The coefficients aj can be rescaled by choosing a different t(r) which can 
be fixed by requiring the conventional parameterization ai + 02 + <^3 = 1 = ct? + «2 + '^3- 
As usual, these equations can be solved only if one coefficient, say ai, is negative and the 
other two are positive. Correspondingly, one direction, the first, contracts whereas the other 
two expand toward larger time. When we approach the classical singularity at t = 0, space 
shrinks only in two directions while the third one expands unboundedly. The total volume 
however continues to approach zero according to 010203 oc t^i+^^a+as = f_ it has been 

shown that the Kasner singularity disappears when the model is quantized along the lines 
of loop quantum cosmology. 
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For Bianchi IX the evolution of the three triad components can be described as motion 
in a non-trivial potential given by 

wip'.p^p') = 2{py(rir2-r3)+pV(rir3-r2)+pV(r2r3-ri)} (15) 

which has infinite walls at small owing to the divergence of the spin connection compo- 
nents. The evolution can then be described approximately as a succession of Kasner epochs 
with intermediate refiections at the potential 2y| . For the Bianchi IX model the refiections 
never stop and the classical evolution is beheved to be chaotic j^]. 

B. Loop Quantization 

Diagonal homogeneous loop quantum cosmology is first formulated in the connection 
representation where an orthonormal basis is given by the p'^-eigenstates 

|mi, m2, m3) := |mi) (g) |m2) C?) 1^3) (16) 



with 

, , ^ exp(kmc) , , 

(^l^)= • (17) 

V2sm(ic) 

Here we are already restricting ourselves to a separable subspace of the kinematical Hilbert 

r~i 

space, which is non-separable |i2S;]. All the information about the classical singularity is 
already contained in wave functions restricted to this subspace such that it suffices for our 
purposes. 

The eigenvalues of the triad operators can be read off from 

p^|mi, m2, m^) = \'^l]>mi\mi, m2, m^) . (18) 

Using the basic operators one can define the volume operator V = y^\p^p'^p'^\ which will 
be used later. Its eigenstates are also \mi, 1712,1713) with eigenvalues 

V{mi,m2,m3) = (|74)^ Vl'^i'^a'^sl ■ (19) 

A kinematical state \s) is described in the triad representation by coefficients Smi,m2,m3 
defined via 

\s) = ^ Smi,m2,m3\mi,m2,m3) . (20) 
mi,ni2,m3 



For a state to be gauge invariant under the residual gauge transformations changing the sign 
of two p-components simultaneously, the coefficients Smi,m2,m3 have to satisfy 

"Smi,r7i2,r7i3 ■5—^1, —m2, ma Smi,—m2,—m3 S—rai,m2,—m3 • (^-^) 

These states are left invariant by the gauge invariant triad operators \p^\ and the orientation 
operator sgn(p^p^p^). In calculations it is often easier to work with non-gauge invariant 
states in intermediate steps and project to gauge invariant ones in the end. 

Together with the basic derivative operators we need multiplication operators which 
usually arise from (point) holonomies hj = exp(c(/)A}ri) = cos(|c/) + 2A}rj sin(ic/) with 
action 

cos(ici)|mi,m2,m3) = |(|mi + 1, m2, ma) + |mi - 1, m2, ms)) (22) 
sin(ici)|mi,m2,m3) = -ii(|mi + 1, m2, ma) - |mi - 1, m2, ms)) (23) 

and correspondingly for C2 and C3. 



1. Inverse triad operators 

From the basic operators we can build more complicated ones. We will later need a 
quantization of the spin connection which is a composite operator containing several triad 
operators. In particular, it also contains inverse powers of triad components which classically 
diverge at the singularity. Since the triad operators have a discrete spectrum containing 
zero, they do not have an inverse. However, general methods of quantum geometry and loop 



quantum cosmology 



29| imply that there exist well-defined operators quantizing inverse 



IpT' = (/7-)^({ci, bY})^ = (/7-)^ ( ^^PM^ih^ ] (24) 



triad components. To obtain these operators one makes use of a classical reformulation, e.g 

:,A\nhi{h^\\p^[ 
ij(j + l)(2j + l) 

which can then be quantized to 

\p%J = ii^ [^ipjij + l)(2j + l)/3]^ {tT,A\TMK\ IpY])^ • (25) 



Here we have indicated that there are quantization ambiguities j30|,|31| when one quantizes 
composite operators. The two most relevant for this paper are indicated by the subscript 
(j, Z). The half integer j appearing as a subscript of the trace corresponds to the choice 
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of representation while writing holonomies as multiplicative operators. The real valued 
/ G (0, 1), labels the various classically equivalent ways of writing \p\~^ as a Poisson bracket 
of positive powers of \p\. 
This operator acts as 

\p^Ji\mi,m2,m3) := {2l)~ {^'jQ ^ A/'/''/j.i(m/)|mi,m2, ms) where, (26) 
AG := ^(i±ip±il and 




The discrete values fj,i{jn) decrease toward lower values for m < 2j j^^. Thus, one can 
see that the classical divergence of the inverse of \p^\ at vanishing is explicitly absent in 
the quantized operator. This will be seen to have further consequences for the approach 
to the classical singularity. Furthermore, the state |mi, 7712,^3) with m/ = (on which 
the classical inverse triad would diverge) is actually annihilated by \p^\'jj due to fj,i{0) = 
independently of the ambiguities j and I in their allowed ranges.^ This allows us to define 
the inverse triad operator (not just its absolute value) by ip^)jj ■= ss^{p^)\p^ji without 
ambiguity in defining the sign of zero. These inverse triad operators will be used later to 
find well-defined operators quantizing the spin connection. Their eigenvalues are 

(/);;/ 1 mi, ma, ms) = (2/)^(i7£j)-Wj~ sgn(mj)/jj(m/)|mi,m2,m3) . (27) 

It turns out to be convenient to define another function Fj^i{m) as: 

F^^iimj) := 2j(2/Ar,)T^/,/m,) , (28) 

so that the eigenvalues are obtained as. 



(p^)^./|mi,m2,m3) = (t^I^I) (2j) ^ sgn(m7) Fjj(m/) |mi, m2, mg) . (29) 



Then, for large j, using also JVj Fj^i{mi) = Fi{ni) actually turns out to be a 

function of /i/ := ^ and with no dependence on j. Explicitly 



30 



3l| 



Fi{q) := (2/(/ + l)(/ + 2)/3)^ [(/ + !) {(g + l)'+2-|g 



^ The usual case with eigenvalues mj^ , which corresponds to a discrete version of the Wheeler-DeWitt 
operator, could be obtained formally from (|26|l with the values I = 2, j = 1/2 which are not allowed 
here, and which would not result in an operator annihilating the states corresponding to the classical 
singularity. 
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-(/ + 2)q{{q + - sgn(g - - (30) 
q-^ (g>l) 

J + 

/,,(m,) ^ (2'+V3)^J^F.(/i/:= t) + o(j^) . (32) 



(0 < g < 1) , and (31) 



There is no suffix j on Fi since it is manifestly independent of j. The hmiting forms 
show that for /i/ ^ 1 one has the classical behavior for the inverse triad components and 
the quantum modifications are manifest for fij < 1 (or for rrij < 2j). 

The inverse triad operators are used in defining the spin connection operators given in 
eq.© (specialized to Bianchi IX) and the corresponding potential (fTHjl . These are diagonal 
in the triad representation with eigenvalues given by, 

r/(/i) = ^ {fiKSgn{fXj)Fi{ijj) + njsgn{iJK)FiilJ'K) - fJ'jfJ'KF^ifJ'i) } (33) 
= 8f [fiifiATjTj - Tk) + cyclic] (34) 

Note that the Tj, as functions of fij, are both dimensionless and independent of j, while the 
potential explicitly scales as j^. For /i/ ^ 1, the F/ become homogeneous functions of degree 
zero of their arguments, an observation which will be useful while discussing the effective 
Hamiltonian. To avoid proliferation of labels, we have suppressed the ambiguity label / on 
the spin connection but have displayed both j, / on the effective potential. 

2. The Hamiltonian constraint 

The quantization of the Hamiltonian constraint is somewhat more complicated in the 
presence of a non-zero spin connection. The details can be found in . Here we summarize 
only the salient points. 

The key difference between the full theory and the homogeneous models with non-zero 
spatial curvature is that the spin connection is a tensor due to the restriction on the co- 
ordinate transformations preserving the homogeneity. Consequently, it cannot be made to 
vanish unlike in the full theory. There is thus no guidance from the full theory regarding 
the incorporation of the spin connection. This loss of guidance is compensated by invoking 



two criteria of admissibility of quantization of the Hamiltonian 
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Recall that the underlying discrete nature of loop quantum gravity implies a difference 
equation form for imposition of the Hamiltonian constraint. The order of this equation de- 
termines the number of independent solutions. Among these, there must be those which 
are slowly varying over the discreteness scale which should correspond to the familiar con- 
tinuum, classical geometry. In essence, the two criteria of admissibility stipulate that the 
difference equation implied by a quantization of the Hamiltonian, must admit a continuum 
approximation which in addition must be (locally) stable in a certain well defined sense. In 
|l9| , a quantization of the Hamiltonian constraint for all diagonal homogeneous models was 
proposed and shown to satisfy the admissibility criteria. 

The quantized Hamiltonian operator is given explicitly as: 

H = 4^(7£j/s:)-i^e^-'^tr[ [r%{A,T)hj{A,T)ChM,'^)r\h{A,T))-' 

UK 

-2(f,fj-n^fi)A,Aj}|/i^(A) [/i^^(A),\/]} ] (35) 

where 



hj{A,T) 



(36) 
(37) 



and the (p^)^^ are defined in equation (j2Zj). Note that f / operators are diagonal in the triad 
basis with eigenvalues denoted by Vi{mi^m2,m^). 

In the presence of non-zero spin connection it is convenient to define a new set of basis 
vectors: 



(38) 



Denoting |s) = Srrai,m2,m3l'^i, '"^2, '"^3) , H\s) = 0, leads to a partial difference equation. 
This partial difference equation is quite complicated and is detailed in jl^ . For the purposes 
of this paper we only need to look at the continuum approximation. 

III. EFFECTIVE CLASSICAL FRAMEWORK 



The effective classical Hamiltonian can be derived from the underlying quantum equation 
in two steps. First one derives a continuum approximation and in the second step derives 
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an effective classical Hamiltonian by considering a WKB ansatz for the continuum wave 
function. 

First, the fundamental difference equation is specialized for a slowly varying discrete 
solution in the continuum regime defined by mj := /{^'jip) ^ 1- The leading terms are 
then independent of the Barbero-Immirzi parameter 7 and lead to the usual Wheeler-DeWitt 
equation for an interpolating continuum wave function 2ll |. 

Explicitly, for a slowly varying discrete wave function, 5mi,m2,m3 := S{p^{mi)), p\mi) := 
^^dpinj, the interpolating continuous wave function S satisfies the Wheeler-DeWitt equation 



p QpiQp2 + cyclic 



+ Wj^i {p^ ,p^,p^)^/ |pip2p3 \ S{p\p'^,p^) 



= - K\pYpTp"''''""{p\p^p)s{p\p\pn- (39) 

Note that in the above equation, the p^ dependence is obtained by replacing m//(2j) = 
A'-i' = P^ / (J7^p)- This will be important in the context of the potential and the matter terms 
(see below). 



In the second step, one takes the slowly varying wave function T := y 5* to be 

of the WKB form T ~ e''^''^. Then the real part of the equation, to the leading order, h^, 
gives a partial differential equation for the phase A, with the phase appearing only through 
its first derivatives with respect to p^ . Eq. (fTUI) suggests the identification p^§^ ■= ^tt/. 
Replacing the derivatives of the phase in terms of tt/ in the partial differential equation 
converts it into an expression which has exactly the same form as the classical Hamiltonian 
f)ll|) . except for the potential carrying the labels j, I. In short, we obtain a Hamilton- Jacobi 
equation for the phase. Classical actions (corresponding to the Hamiltonian) along dynam- 
ical trajectories are well known to solve the Hamilton- Jacobi equation. Thus constructing 
slowly varying solutions of the Wheeler-DeWitt equation naturally leads to analyzing the 
trajectories determined by the classical Hamiltonian derived in the above manner. 

The classicality of the Hamiltonian presumes that in (jH^ the potential and the matter 
term are of order hP. This feature however depends on the values of p^ /{j'jip), because the 
ambiguity labels j, I control the p^ dependence of the potential. For instance, taking the 



In the equation (65) of this reference, the coefficient of the first terms should be 1 instead of j 
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value of j to be large, allows us to divide the mj ^ 1 continuum region into (i) j ^ mj ^ 1 
and (ii) mi ^ j ^ 1 sub-regions. Since Fi is a function of m/{2j), the two sub-regions are 
more conveniently demarcated as <^ J7^p and ^ J7^p respectively. 

In region (ii) the eigenvalues of the inverse triad operators correspond to the classical 
values (inverses of triad eigenvalues) and consequently the spin connection acquires no de- 
pendence on 'jip when /// are expressed in terms of the p^ . The potential then reduces to 
the exact classical expressions without any dependence on •yip. 

By contrast, in region (i) the Fi{fij) have a power law dependence with positive power. 
This makes the spin connection have a positive degree dependence on its arguments and 
introduces inverse powers of 'yip in the potential (and matter) term when fij are replaced 
by p^ /{jyip)- These terms clearly reflect non-perturhative quantum effects in region (i). 

Thus the argument that the effective Hamiltonian derived via the WKB ansatz is classical 
(hP) is strictly valid only in region (ii). Noting that in region (i), the inverse triads go 
to zero as p^ 0, we will assume the validity of the effective Hamiltonian in region (i) 
even though the potential has inverse powers of yip. In this manner, in an approximate 
sense, we incorporate non-perturhative modifications of the potential in an otherwise classical 
description. 

The validity of the classical description is controlled by the validity of the WKB approx- 
imation while how much of quantum modifications can be seen in the classical framework 
is controlled by j. Similar remarks apply to the matter Hamiltonian terms since the matter 
density also typically gets quantum modifications from the inverse volume operator. 

The classical Hamiltonian deduced from the WKB approximation in the above manner, 
including the extension to region (i), is what we will now refer to as effective Hamiltonian. 
Notice that imposition of the underlying quantum Hamiltonian constraint has automatically 
implied that the effective Hamiltonian must be zero which we interpret as the (modified) 
classical Hamiltonian constraint. 

To study the dynamical implications of the effective Hamiltonian constraint, recall first 
that in general relativity we can have two distinct notions of dynamics or time evolution (see, 
e.g., One is evolution with respect to an external time which is a suitable coordinate 

in the space-time description of the geometry and the other is evolution with respect to 
an internal time which is one of the degrees of freedom treated as a 'clock'. The latter is 
obtained by solving the Hamiltonian constraint for one suitably chosen momentum variable, 
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say some ttq, which becomes a function of the remaining phase space coordinates and the 
conjugate coordinate g° identified as the internal time. The evolution of the remaining 
degrees of freedom with respect to g° is specified by the new Hamiltonian, H := — ttq, 
which is in general time dependent. Given a solution of the internal dynamics, one can 
systematically construct a corresponding solution of an external dynamics. What a quantum 
Hamiltonian constraint does is to provide a modified classical Hamiltonian constraint - our 
effective Hamiltonian - via the continuum approximation as outlined above. This results 
in the modification of the trajectories of the internal dynamics and correspondingly of the 
trajectories of the external dynamics. In the context of homogeneous models, the internal 
dynamics turns out to be time dependent but an un- constrained system. 

In the subsequent sections, we will analyze the internal dynamics determined from the 
modified Hamiltonian constraint. 



IV. INTERNAL TIME EVOLUTION 



We begin with the effective Hamiltonian (without matter) obtained from the continuum 
approximation, 

KNHfi{-Ki,q') = -l(7ri7r2 + 7r27r3 + 7r37ri) + l^,,Ke2«\e2<?',e2'?') , {g^ vr^} = (40) 

Here, we are using dimensionless variables (or equivalently set the unit of length such 
that = 2) and have put raj := e^^^ > 0. The potential is given by eq. ()34|) with the spin 
connections given in eq. (j33|) . 

First we diagonalise the kinetic term in the Hamiltonian. 



[7ri7r2 + 7r27r3 + ttstti) := 



1 



A :-- 



( 



1 1^ 

1 1 
1 1 Oy 



^ — 

-1 1 1 
\/6 

. ^ L J_ , 

\ V6 ^2^3/ 



4 

1 \ / 



-t{^ ODO^ where, 
4 



(41) 



v 



-1 





o\ 


/ 2 


1 


1 \ 





-1 








1 


1 










V 73 


1 


^ / 



(42) 



Defining tx' := O'^n, q' := O^q we get the kinetic term as — ^tTq + |(vr^ + tt^) and the 
potential is obtained as a function of q' =: (g+, q^ , g°) by substituting q = Oq'. Explicitly, 



2 



1 J-gO J_gi 



2j 
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Clearly, g° = ;^(?^ + + Q^) — In (volume). For an internal time description, this is 
a suitable variable. The internal time description is obtained by solving if|f = for ttq: 



-no^n{q'^,q ,q ,7r+,7r_) := 



(44) 



H. is to be used as the new Hamiltonian for analyzing the behavior of as the volume 
shrinks to zero monotonically i.e. as (f ^ — oo. Note that there are no longer constraints 
but the new Hamiltonian is explicitly 'time' dependent. We are also suppressing the labels 
J, I on Ti and h. 

The equations of motion then take the form, 

dg± 1 d7r± 1 dWu 



^± ' = ^ (45) 

(46) 



2^ ^ ' d?o ^ dq^ 

dV _ 1 dWj^i IdWj^idq^ 



d(g0)2 2h dq^ h dqO dq^ 
Along a solution of the equation of motion, we can eliminate h using h — 2h{{q'^y + 
(q-f) + 2Wj^i or /i = 2Wj^i/(l - 2((q+)^ + (r)^))- This leads to 

dV_ /,_2f^>;y_2ftVU 1 ^+ 1 ^^4^ (47) 



d(g0)2 I \dqy VdgV I {"^Wj^i dq^ 2Wj,i dq^ dq^ 

It will also turn out to be convenient to introduce polar coordinates in the q"^ anisotropy 
plane. This is easily done: Define q'^ :— Rcos{9), q" :— Rsm{9). The conjugate mo- 
menta then are given by 7r+ := 7rijcos(^^) — ng^^^^, 7r_ := 7rijsin(6') + ng ""^^^ . The squared 
Hamiltonian, h, and the equations of motion are given by 

h{R, e, nn, no, g°) = ^ K + R-'^l) + 2Wj,i(p\q', R, 9)) , (48) 
H(i?,^,7ri^,7r,,g°) = ^yh{R,e,7rR,7r0,q°) . (49) 

Denoting the d/dq^ by an over dot, the Hamilton's equations of motion take the form, 

(50) 
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The second order equations are obtained as, 

Along the dynamical trajectories, we can eliminate the conjugate momenta from the 
Hamiltonian and obtain, 

1 - 2 + RHA 

h^^ = > (strict inequality) . (54) 

2ryjj 

When the potential is zero, as in Bianchi I, the right hand sides of fl52l I53|) are zero. It is 
easy to see that the Kasner solution is obtained a.s 9 = constant, R{q^) = + constant. 

When the potential is non-zero, the equations are quite complicated. To proceed with the 
qualitative analysis, we will study the effective potential in the next section. 

V. STUDY OF THE EFFECTIVE POTENTIAL 

The effective potential is more conveniently studied in terms of the polar variables. The 
following expressions are useful for our subsequent analysis. 

:= R cos 9 q~ := R sin 9 

^ -A= R COS 6 + R cos 9 

fii = Qe^ = Qe ^ 

/i2 = tt/e ^ = Qe ^ ^ ^' (56) 

^ --2=i?,cos6»(l+V3tan6») ^ - ^i?,cos(6»- f ) 

/^s = ^ = Qe ^ ^ 

Clearly, /ii/i2Ai3 = is independent of the g^. 

Using the obvious abbreviation, Fj := Fi{ni), (label / on F is suppressed) the potential 
is given explicitly by. 



2f 



[-2 {/ilAi2(Ail^2 + fJ'2Fl) + /i2Ai3(Ai2^3 + /i3^2) + Ai3Ail(/i3^1 + /^l^s)} + 

3 {{fiifi2F3f + {fi2fi3Fif + {fi3fiiF2f} + {/iiFaFs + fi2F3Fi + fi^FiF^} 
-Q' {f^ifi2F3{F^ + Fi) + fi2fi3Fi{F^ + Fi) + fi3fi,F2{Fi + F^)} 

{fiifi2Fi + fi2fi3F,' + /X3/iii^|} + Q' {F^F^ + F^Fl + F^Fl]] (57) 
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We define the far regime, by /i/ ^ 1 and the near regime by < ^ 1. Note that the 
regimes are delineated in terms of the fij and not in terms of the anisotropy variables q^. 
The regimes are realized for any bounded region in the anisotropy plane {q^- q^ plane) and 
sufficiently large volume {Q ^ 1) and sufficiently small volume {Q <^ 1) respectively. In 
these regimes, the potential simplifies to the forms, 



/^3 



2-1 
l-l 



2-1 

759+(t^) cosh I T-^V2q' 

V2q 



2-1 
. 1-i 



-my 



2-1 
' l-l 



l+l 



1-' J^„o,±z2L\ 



^ ^ ^ cosh 



l-l 
2-1 



l-l 

+ e^^^^^^ cosh f 2v^g" 



(58) 



fJ'3 



1 -^gO 



e 



- 4e"^''^ cosh (^2V2g"j 

+ 2eV5^^ I cosh [AV2q-^ - l| j 



(59) 



Comparison of expression (jSHl) with Misner's expressions j26[ suggests the identification: 
g+ , q- , g° - VSfi. We will continue to work with the g's. 

Note that although in the two limiting cases the g° dependence factors out in the potential, 
this is not true in general. This is different from the classical potential which is just the far 
form of the effective potential. 

The factorization property of the far and near potential allows to show that there can 
be no bounce in the corresponding regimes, i.e., once the volume decreases, it must keep on 
decreasing until it reaches zero [2^ . This follows from the external time dynamics generated 
by the constraint H = — 7rQ + |(7r^ + 7r^ ) + 2Vr = 0. For g° and its conjugate the Hamiltonian 
equations of motion are g° = — 27ro and ttq = —dH/dq^ = —2dW/dq^ = —2kW with positive 
constants k in the near and far regime. This implies 



dt 



-kq°/2 



ke-'^^'/'inl - 2W) = ke-^^^'^Til + 7r!)/2 > . 



Thus, once the volume decreases, i.e. g° < and ttq > 0, it cannot turn around since hqC ^'fl'^ 
starts out positive and cannot decrease. This proof, however, relies on the factorization of 
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the potential into a part depending exponentially on g° and a part depending only on the 
anisotropies. In regions where such a factorization is not possible, it is not clear if bounces 
or oscillations are ruled out. The modification of the potential could lead to a bounce at 
positive volume as observed in for a closed isotropic model with scalar matter. If a 
similar bounce would occur in the effective Bianchi IX dynamics, the system would clearly 
not show the classical chaotic behavior. Thus, we can assume here that there is no bounce 
in order to investigate the potentially chaotic behavior. 

Other asymptotic regions correspond to two of the /i/'s large (small) and the third one 
small (large). These are realized for example, for a fixed and moderate Q and large R. The 
behavior of /i/'s is then controlled by the angles. These regions divide the anisotropy plane 
into angular sectors with R large. In these regions, we could use the asymptotic forms of 
Fi{fii)^s and obtain the corresponding approximate forms of the potential. These forms are 
of course different from the near and far forms displayed in equations ()58l 13^ . These are 
discussed in more details in section IVII 



A. Approximate wall picture 

One can see the qualitative behavior of the potential (classical or effective) by specializing 
to = or 6^ = Oandvr. For this q~ , we get U2 = 113 := a and Ui := u = %. It follows 
that r2 = T3 and the potential simplifies a bit. A wall, i.e. positive potential, is seen for 
/i ^ 1 so that Fi{fi) ^ fi~^ may be used but u is not so small that Filu) ^ 4i/^ can be used. 
Making this substitution for -Fz(/i) one obtains, 

= ^ i^FKz/) , Ti ^ 1 - ^fi'F^iu) ; 
iy,-z(z/,/i,/x) ^ 2ffi'F^{u){3-2uFi{u)}-8ffi' 

^ 2f%F^{u) {3 - 2uF,{u)} - 8f^ , (60) 

= : Vi{Q,u)- 8f-^ (definition of Vi{Q, u)) . (61) 

We refer to Wj^ii^u, ^, ^) in eq.flUUj) as the wall potential and to Vi^Q,^) defined in as 
the approximate wall potential. 

The last term in the wall potential can be neglected for u ~ o(l) and Q^/iy = fi^ ^ 1, 
since F;^(z/)(3 — 2z/F;(z/)) is of the order of 1 for this range of values. Neglecting this term and 
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substituting = (volume)'*/(2j)^ in the approximate wall potential, gives the expression 



quoted in |22|. 

For z/ ^ 1, one can replace Filu) by u^^ as well leading to VilQ,!/) ^ Ij^Q^jv^ ^ 
^ (volume) ^''^e~vl'' matching with the leading term of eq. (j59j) for g~ = which is the 
classical wall potential. 

In contrast to the classical wall, the wall of the effective potential has finite height. The 
location of its maximum along the g+ axis is defined by = 0. Its solution will define 
a function q'^{Q). The motion of this position is obtained by solving ^(^) = 0. Recall 

4 + 

V = Qe^'' =: Qe^' such that 

dVi{Q,u{Q,x)) _ ^ ^ ^ = 
dx dv 

d^Vi d^Vi , , 

Thus, a (local) maximum of the approximate wall potential can be defined either in terms 
of derivatives with respect to a; = -^(l^ or with respect to v. 

Notice that for the classical wall, there is no local extremum and therefore one has to 
track the motion of a fixed value of the potential. Thus, when the volume decreases, 
must also decrease., i.e. the wall moves away from the isotropy neighborhood (note that the 
wall is seen for /i ^ 1 which implies < 0). By contrast, the approximate wall potential 
does have a local maximum at some i^o- Fixed v implies that, with decreasing volume, the 
height oi the approximate wall decreases and q^ increases. Thus the approximate wall moves 
towards the isotropy neighborhood with decreasing height. Eventually of course the wall 
potential must be used and the height actually crosses zero for a certain volume depending 
upon the values of j, /. 

It is difficult to obtain a convenient approximation for a wall along an arbitrary direction 
and numerics must be used. Numerically, the full effective potential becomes zero every- 
where, in the anisotropy plane, for Q ^ 1.08511 or for the dimensionless volume of about 
(2.172j)3/2 (for the choice / = i). 

From these studies, it is apparent that the principal reason for the usual chaotic behavior 
disappears suggesting that the usual 'Kasner epochs' will continue without any replacements 
once the volume becomes comparable to the Planck volume. This in turn suggests that the 
Kasner solution should be stable under the effective potential. 
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B. Stability of Kasner solution 



This is most easily done in the Cartesian form (jUj). Substitute q'^iq^) = -ft'^(g°) + (5g^(g°) 
where := ~^u^q^i '■= cos(0),m~ := sin(0). Linearization gives, 





AWj^i dq^ ^ 2Wj,i dq' 
1 dWj^i 1 dWj 



Kasner 



^W,,i dq^ 2V2 Wj,i dq^ 



Kasner 

(63) 
(64) 



It is obvious now that the determinant of the matrix is zero and the only (possibly) non-zero 
eigenvalue is given by the trace: A := k+ cos0 + k~ sin0. 

The eigenvector corresponding to the zero eigenvalue is given by 5q^ = — cot0(5g+. The 
first braces in the middle equation of (jUHj) then vanishes. This zero mode implies that the 
perturbed solution is also a Kasner solution. But we know that the dynamical equations do 
not admit Kasner solution since the potential is generally non-zero and h is to be finite and 
positive and this solution of the perturbed equation must be rejected. 

Thus the Kasner solution is stable if A > (positive because we are considering the 
backward evolution, decreasing q^ ). Up to here, no specific form of potential has been used. 

To compute we need to compute the potential on the Kasner solution. Substituting 
Kasner K^{q^) in the gives, 

2j^, = e^.°(l+|-s^-^sin^) ^g5^ 

It is obvious now that all the exponents are non-negative over all implying that all fij will 
(generically) decrease as g° is decreased. Only for (p = 120° where fi2 = = 240° 

when /i3 = l/(2j) and = 0° where fii = l/(2j), all exponents are positive. By taking 
j > 1/2 we can ensure that all the /i/ are much less than 1 at some initial g°. We can then 
use the near-form of the potential to compute the derivatives. 
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Defining, generally, 2j/i7 =: e'^'^'^^'^ ^ one can obtain the derivatives of the potential as, 

= ^ (66) 

= -if (jI^) " Wii^i + A) + i4(.i4 + A-?) + + i^l)] (68) 

where L = The logarithmic derivatives of the near potential are given by, 

1 dWj'P'' _ 2(L + 2) ^ 2-1 

W]f^~d(f~ ~ 1 - / 

1 5iy°f^'' 2 -1 

+ 2)/i^ + (4 - + /i^((4 - L)f^l - (L + 2)/i^)] (70) 



(69) 



1 c^l/l/"'''^'' 



+/^^((L - 2)/i^ + L^i) - ^^(L^f + (L - 2)/x^)] (71) 

These now have to be evaluated on the Kasner solution. 

Observe that the logarithmic derivatives of the potential with respect to are ratios of 
homogeneous polynomials of the same degree. Hence any common factors in fij will cancel 
out. From it is obvious that effectively we can take, 

/i, := ^-3, := u-\ /is := where /i := e^'°'^ u := e''"^''^^ (72) 

The common factor canceled is Q{q'^,j)fJ'- Substituting these in (fflU l7T |l leads to. 



1 fjW^'^^''^ 9 



[2(L - l)/i-=^^(i/' + z/-') + (4 - L)/^-6(z/^ + z/-^)- 

(L + 2)(z/^-2^z/-(^-2) 

=: ^^~'«+ (73) 
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Ji'- 



linear Q - 



= !(£±l)e„3^-^^;± (75) 

2(^^ + 2) 1 a_ 

2(L + 2)-^±^ + ^^l (77) 



A 



1 

7i 



D D 

The eigenvalue A is a function of and gf°. Generic stability of Kasner solution requires 
A to be positive for all 4> and at every less than the initial q^. A necessary condition 
for this to be true is that a\ + < 4{L + 2)'^D^. This is easily checked independently of 
(j), completing the proof that the Kasner solution is stable under inclusion of the effective 
potential for all (f) and for all g° less than the initially chosen g°. 

An identical calculation can be repeated for the classical potential by using the far form 
of the potential (which is independent of Z, of course). However, now A does not have a 
definite sign for all (f>,q^. For a fixed 4>, A changes sign. The positive sign corresponds to 
'Kasner epochs' while the negative sign corresponds to 'replacement of Kasner epochs'. 



C. Relation to Kasner Sphere 

It is conventional to analyze the asymptotic behavior of the approach to a singularity in 
terms of the scale factors of the classical geometry. We can obtain such a description by 
noting that (all scale factors, triad components positive), 

= ajGK (78) 

volume /^.viQ^^^ . 

ai ^ — = 2j^^ I.e. 79 

— 7^ -Rcos 9 /orv\ 

ai — e ^ (80) 

a2 = g^g+^i?cose(l-V3tane) ^g^^ 

For the Kasner solution, -R(g°) = "71^°'^ ~ ^° while the standard parameterization is 
ai ~ X^'^^ = g2a,int g^pjj i-j^j^^ ^^aj = 1 = XI/ «7 ('Kasucr sphere'). Putting 5° = 2Clnt, 
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comparison implies, 

CKi = {1 + 2 cos ^0} 
v3 

= |l — cos 9o{l — Vst&n Oq)^ 
0:3 = {-^ ~ + v^tan ^0)} (83) 

The Kasner Sphere conditions are then automatically satisfied with the choice ( — ^ 
which identifies the co-moving time t as e^^° which is consistent with the definition of q^. 

It is obvious now that a point on the Kasner sphere corresponds to a choice of 9q. Oscilla- 
tory behavior corresponds to a motion of the system point passing through different angular 
sectors mentioned above. 

Observe that a deviation from Kasner solution {9 non-constant) does not imply oscillatory 
behavior. The 9 motion may still be confined to one of the angular sectors leading to a 
Kasner-like behavior. For the classical potential with its confining nature (infinite walls), 
non-constant 9 implies traversing all the angular sectors and thus oscillatory behavior. For 
the modified potential such an implication is not automatic although conceivable. 

In the next section on the qualitative analysis we will see that the asymptotic solution is 
indeed a Kasner solution. 



D. Domain of Validity 

Before we consider the qualitative analysis, it is useful to make explicit the domain of 
validity of the effective Hamiltonian. Recall that the derivation of the effective Hamiltonian 
is based on the validity of the continuum approximation and the validity of the WKB 
approximation. 

The continuum approximation itself is valid for m/ ^ 1 which implies /ij ^ e := l/(2_7). 
In terms of the q^'^ these restrictions translate to: 

Q > e <^ g° > ; 
^ ln(e/g) « q+ « ^ ln(g/e) ; (84) 



]r\ O A — n+ \ <!<r n- ^ 



ln(Q/e) -q+\ < g- < ^ ^ln(Q/e) - q 



2 
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Equivalently, using Qle = Vol^/^ we have, 



1 




(85) 



(86) 



For a given /, the effective potential becomes negative for some Q^, independently of j and 

3 

the corresponding vol^ ~ (2j(5^)2. To be consistent with the domain of validity, Vol^ > 1 
must hold which can be arranged by choosing j sufficiently large. Having chosen Vol;*, the 
limits on follow. 

Some typical values are: Qq.i ~ 0-91 ! Qo.s ~ 1-086 ; Qo.9 ~ Clearly, for smaller 

/ we need to choose a larger j in order to stay within the domain of validity. (In this 
context it is interesting to note that both conceptual considerations in the full theory and 
)henomenological studies in the context of early universe inflation point to values j > 1/2 
Thus we see that disappearance of the walls can be se en g enericallv and reliably within 
the effective Hamiltonian picture. For the plots shown in j2^, the ranges of are in the 
domain of validity. 

For the stability analysis above, for large enough j, we can have all /x/ -C 1 even though 
g° > justifying the use of the near form of the potential. The Kasner stability is thus seen 
staying consistent with the validity of the effective Hamiltonian. 

One reason for this demarcation is that for any given fixed volume, the large R form of 
the effective potential is very different from that of the classical potential. In particular, the 
classical potential is always large and positive while the effective one is large and negative. 
Only the inner side of the walls of the effective potential 'matches' with the classical potential 
for large volume. Consequently, the 'initial conditions' for the internal dynamics can be 
either in the inner side of the walls or the outer side of the walls - an option not available 
for the classical potential. 



In the next section, we will go beyond the domain of validity by considering the limit 
Vol ^ 0. 
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VI. QUALITATIVE ANALYSIS 



In this section we will pretend that the effective description is valid all the way to zero 
volume or g° — oo. Although we know that this cannot be true, we are considering this 
by way of contrasting the dynamics implied by the classical potential. We will therefore 
be interested in the asymptotic behavior as g° — oo of the trajectories in the anisotropy 
plane. This analysis proceeds in stages. Asymptotically one of the following two possibilities 
can happen: either R{q^) remains bounded or it does not. In the first stage, we will first rule 
out a bounded asymptotic motion while in the second stage we will show that asymptotic 
motion always approaches a Kasner trajectory. 

Proposition 1 Asymptotically, R{q^) cannot remain bounded. 

The proof is by contradiction. Assume that the asymptotic motion is bounded. Then all 
the /i/ must become vanishingly small due to the explicit exponential g° dependence. Thus 
for sufficiently early time (small volume), we can use the near form of potential valid for 
/i/ < 1. 

Consider now just the radial motion, R{q^). Since it is bounded, it can either tend to a 
limiting value R^, or it will keep oscillating. If it has a limiting value, then we must have 
R, R = at R = R^:. However, for the near form of the potential, its radial derivative is 
negative. The radial equation of motion ()52|) then implies a contradiction unless possibly 
R, = 0. 

For i? ~ 0, we have, 

1 

P^?r (i?, 9) ^ -2(2jr ^ e75(^)'?° (s + ^^f^R') . (87) 

Since at i? = the potential is negative and must be strictly positive, we must have 
(i?^ + R^O"^) > 1/2. But this is impossible at R^ = 0. Thus, R{q^) cannot have any finite 
limiting value. 

If R{q^) oscillates, then it must have several extrema. Now the radial equation (|52p 
together with the fact that < implies that R must be positive at any extremum. 

Thus the radial motion can have at the most one minimum and oscillatory motion is ruled 
out. Since both the oscillatory and the limiting radial motions are ruled out, we conclude 
that the trajectory cannot be confined to a bounded region and the proof is completed. 
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For analyzing the case of unbounded motion, it is convenient to use a vector notation: 

^ = _L =: V I = 1, 2, 3 





1 


Pi 










8 




3 



I 

In addition, we also have the expression for fij in terms of R, 9 given in equation ()56p . 

Observe that the ratios of the /i/'s do not have any explicit dependence on g° and since 
R{q^) is very large, these ratios either diverge, vanish or become equal to 1 (the last case 
only along three special directions). Values of individual /i/'s however depend on the relative 
growth of R{q^) and Ig*^]. We will say: R{q^) is dominant if -^j^p ~^ oo, sub-dominant if 
— i> and marginal if -^j^ where ^ > is constant such that ^ = 1 corresponds 

to a Kasner motion. For sub-dominant R{q^) and marginal one with ^ < 1, all /i/'s vanish 
(generically in the marginal case) while for dominant and ^ > 1 cases, two of the /i/'s diverge 
(vanish) while the third one vanishes (diverges) since the product of all /i/'s vanishes. Since 
2j > 1, all /i/ are either much smaller or larger than 1 so that the corresponding forms of 
the F{fij) can be used to approximate the effective potential. To get explicit forms, we have 
to consider various cases. Noting that the potential is invariant under cyclic permutation of 
the /i/, it suffices to consider the cases: (/) /^s ^ /i2 ^ Pi and (/') ij,2 ^ Pi ^ Pi ■ It is 
easy to see that (/) corresponds to large R and a^, ■ q < a-z ■ q < (fi ■ q from which it follows 
that q must be confined to < 6* < |. Likewise, the region (/') corresponds to large R and 
-f < ^ < (see figure □). 

Within these regions (ordering of the /i/'s), we have further sub-cases, 

(a) /i3 < /i2 < /il < 1 , (b) /i3 < /i2 < 1 < /il , (c) /is < 1 < /i2 < /^i • (89) 

The case where all 1 -C /i/, would imply that their product is larger than 1 which is not 
possible for g° — > — oo. 

For these cases, we will derive the asymptotic form of the potential and analyze the 
equations to deduce that only marginal R{q^) is possible with ^ = 1. 

To derive the asymptotic potential, we prove a couple of lemmas: 



Lemma 1 If R{q^) — > oo as g° — — oo then for any real numbers l,m,n, {^i)\^j)'^{^k] 
either diverges, vanishes or goes to (2j)^('+™+"). 

28 



n 




FIG. 1: The asymptotic regions I and I' of the anisotropy plane. R is outside the circle, da varies 
with the ambiguity parameter I staying within I. The vectors shown are unit vectors. 



The proof is obvious since jjii 
either positive, negative or zero. 



~ e^'^ ^^'^ )cos6>7(g ) asymptotically, the exponent is 



Lemma 2 The spin connection components obey the inequalities: 



(a) 1 > /// > //J > //K 
(6) /// > 1 > //J > //i^ 

(c) /X/ > /Xj > 1 > /X/^ 



T/^ > > , Vk:> TjVi ■ 

> Tj > r, , > r^r, . 



To prove the above lemma, it should be noted that the ratio of any two spin connections, 
say Vi and Vj can be written as 

(90) 



rv 



fJ-J 



UK IJ.I 



where the term inside the parenthesis is explicitly in reverse order. Also it can be shown that 
for all cases and for all allowed values of the term inside the bracket asymptotically either 
vanishes or goes to a constant or diverges but the divergence is weaker than the convergence 
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It follows that the potential takes the form Wj^i{fl) ~ Sj'^fiifijTK- From Lemma d it 
follows that Tk will become a monomial in the /x's. The potential then takes the form: 

Wffifl) ^ - Af ( j^] (91) 



/ + 1 



or -4jVj' (92) 
or -4j2|l+(^)T^(2jr^}/^/, 
Wf}{il) ^ - Af^ij' . (93) 

Thus we have proved, 

Proposition 2 In all cases where the {ij either vanish or diverge, the potential necessarily 
takes the form, 

Wj^i{R,e,q^) ~ _Ce^9»+-a (94) 

where C, p, a are constants with C > 0, p > and < 2p^. 

Thus in all the cases, the potential has the form (jMj) . In fact, for the sub ) and 

the first of the sub-case (b) above, we have a^/(2p^) = (3/^ — 6/ + 4)/(3/ — 4)^, which is 

smaller than one for all < / < 1 (and exactly one for the limiting case / = 1). In the 

remaining two of the sub-case (b) and the sub-case (c), we have (?^/(2p^) = 1. Notice that 

■J 

the asymptotic potential is negative. Together with h > 0, this implies that 1 — 2g < 0. 

Using the asymptotic form of the potential, the Cartesian form ()47|) of the equations 
becomes {v := q), 

^'-0(^ + P^) ' "^'^l ' 2p'-^'>0. (95) 

Notice that the above equation is a first order, autonomous equation for the velocity. 
The (vector) equation has fixed points at -y^ = ^ and at v = —j^- For the latter, isolated 
fixed point, "^^^ = |^ < | and hence this is not accessible. 

Taking scalar product of the equation (j95|l with 2v one gets, 

j^^(^^-i)(2,.^)(|-|+l) . (96) 

The inequalities on the iP and imply that the third factor is always positive and 

hence must decrease as g° is decreased. Clearly, iJ^ is bounded above by (wmitiai)^ and 

30 



bounded below by |. Since the equation for is monotonic, asymptotically 'tp must reach 
the lower bound. 

Taking the 'cross product' of the equation with a one obtains, 

This equation shows that either a x t? = or, if non-zero, cannot change sign. If (j) denotes 
the angle between v and a, then the cross product is f crsin0 {v := \v\,a := \cx\) which is 
clearly bounded. Hence it must also approach a limiting value. Since both the magnitude of 
the velocity as well as its component orthogonal to a approach a limiting value, we conclude 
that asymptotically v must approach some {7* such that = ^. This also means that Riq'^) 
must be marginal with ^ = 1, both the dominant and the sub-dominant as well as marginal 
with ^ 7^ 1 possibilities for R{q^) being ruled out. Since in the sub-cases listed in equation 
()89|). (b) and (c) have at least one of pi, fi2 to be larger than 1 and for ^ = 1 marginal 
behavior all /x/'s must be less than 1, these two sub-cases are ruled out. Thus, sub-case (a) 
is the only relevant case and, corresponding to region (I), we have a = ffa which is read off 
from equation (PT|) as, 

aa = ^<?i + 2c?2 . (98) 

As shown in the figure ^ lies within the region (I). 

To study the nature of a candidate fixed point v^: of equation ()95p. linearize the equation 
by putting -u = + 6v. This leads to, 

Sl) = {v^ ■ 6v)A , A := (T + 2pC . (99) 

Clearly, 6v orthogonal to v^: is a zero mode of the linearized equation which does not 
evolve. Taking scalar product with gives, using = |, 

(C ■ 6v) = • 5v){v^: ■ a + p) . (100) 

Linear stability of {Ik thus demands p+cr(cos0*)/-\/2 > 0. Since -\/2p/cr > 1, this condition 
is satisfied for all (p^,. Thus we have proved that the asymptotic velocity under the evolution 
equation (195^ is a Kasner velocity and every Kasner velocity is admissible. Furthermore, 
this velocity is approached exponentially as implied by the stability analysis. 
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One can also derive a relation between an initial velocity and the limiting velocity. Denote: 
f|l := f COS0, v± := t>sin0. Taking scalar and cross product of equation with a gives, 

^1! = (^'-0 + ' = (^'-0(P^^) • (101) 

If fx = initially, then it is identically zero and therefore 0^, = or tt depending upon 
the sign of v\\. The v\\ equation can be easily integrated but we do not need the explicit 
solution. If f _L 7^ initially then it must retain its sign. In this case taking the ratio v\\/v± 
gives a differential equation for v\\ as a function of v± which can be immediately integrated 
to give an orhit equation, 

= A (^v\\ + (102) 

where A is a constant of integration. 

The orbit is a straight line in the velocity plane always passing through the point v\\ = 
— ^ ( > — ) ,'v± = and with a slope determined by A. This point is exactly the isolated 
fixed point of the equation. 

Note that since v± neither A nor {v\\ + ^) can be zero for any finite We also 
know that asymptotically v±_ — > t>_L* := ^sin0* and v\\ := ^cos0*. As noted 

earlier, only is consistent with asymptotic Kasner motion and for this case 2p^ — is 
strictly positive. This implies that (fy* + ^) is non-zero and therefore f_L* is also non-zero. 
Evaluating equation ()102|1 for initial values of {v\\,v±) determines A as a function of v,(f) 
while evaluating the same equation for asymptotic values (f fj,*) determines the angle 0*. 
Explicitly, in terms of A, 



(7 / .o O"^ 



The two roots (± signs) for sin((/)*) always have opposite signs and since the sin(0) does 
not change sign, the initial condition sets the sign and selects a unique limiting value 0=,,. 

Now the trajectory in the anisotropy plane can be obtained trivially as, q{q^) = v^^q^ + q 
which leads to, 

R{q') ~ --^g° + 0(l) , %°) ~ ^o + 0|^-^| , (104) 
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where 60 = 0* plus the angle a makes with the g"*" axis. Notice that asymptotically q points 
in the direction opposite to that of v^,. 

Let us now compare a similar analysis for the classical potential. In the classical case, one 
does not need to distinguish the sub-cases of (jHI?|) as the potential is directly given as a sum 
of monomials in fij. The asymptotic form in region (I) is the same as given in Proposition 
|2]with C < 0, p = :^,<? = —4(73 = (:^'4a/2). This implies two differences compared to the 
effective case: = | < 1 and < ^. The first order, autonomous equation for the velocity 
remains the same as (except for the inequalities) and so do the equations (jHSl EZI) • The 
isolated fixed point of gets ruled out as before. Unlike the effective potential case, the 
third factor in (|96|) no longer has a fixed sign. If it is zero, then f ^ gets frozen at a fixed 
value and if it is non-zero, it cannot change its sign. Since the evolution of f ^ is bounded, 
once again we get that — > although we cannot conclude at this stage that = |. As 
before, v sin is bounded and evolves monotonically thereby reaching some limiting value. 
Thus we conclude that -u -u* as before though without the condition = |. 

Properties of a candidate limiting velocity are analyzed by linearization, v = + 6v. 
This leads to, 

6v = 7A -|- 'fpSv + {v^, ■ 6v)A , A := a + 2pv^ , 7 := . (105) 

Clearly, 5v ^ Q must hold and therefore either 7 = or A = 0. A = is ruled out since 
this gives v1 = 2 and hence 7 = must hold. We have recovered the = 1. The rest of the 
stabihty analysis now gives, exactly as before (p -|- " ) > or cos(0=k) > —\. Unlike the 
effective potential case, now we have a restriction on the asymptotic angles 0*. The angle 
0* can again be obtained as in eq. ()103|1 . however due to ^ = 4, the square roots are real 
only ii <\. 

The orbit equation provided a very simple pictorial representation shown in Figure El 
The two differences between the effective potential and the classical potential control the 
location of the isolated fixed point {v_[_ = 0) and the setting of initial conditions. For both 
potentials the continuum of fixed points is determined by the circle: v"^ = ^. For the effective 
potential, the isolated fixed point lies inside the circle while the initial condition must be set 
outside the circle. For the classical case it is exactly the opposite. The monotonic decrease 
of f ^ in the case of the effective potential and monotonic increase of f y in the case of the 
classical potential shows the direction of traversal of the orbits. The exponential approach 
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to C follows from the linear stability analysis. 




p 




,-1/2 



,-1/2 



A = - 3 



(a) 



(b) 



FIG. 2: The point P and the circle constitute the fixed points. The heavy, directed straight hne 
segments are typical orbits of the velocity dynamics in the velocity plane. The U|| , v±_ are relative 
to a. The angle shown is the total change in the direction of the velocity vector. The figure (a) 
corresponds to the effective potential while (b) corresponds to the classical potential. 

How does one understand the asymptotic Kasner solutions for the classical potential given 
the well known oscillatory behavior? Recall that the asymptotic direction reached in the 
anisotropy plane is opposite to that of ■u* (since g° is negative). The stable Kasner velocity 
directions are confined to — ^ to ^ relative to a which makes an angle | to the axis. The 
corresponding trajectories therefore satisfy < 6q < 2-n which precisely excludes the region 
(I)! The qualitative analysis rests on the assumption that remains valid throughout the 
evolution. This assumption must fail in the case of the classical potential — a trajectory 
that enters an asymptotic region (I, I' or their 120° rotations), must exit it in a finite g° 
interval. By contrast, for the effective potential, since all 0=,, are stable, there always exist 
directions such that Oq is in the asymptotic region. The remaining directions will again 
correspond to exiting the region only to eventually enter another region with compatible 9^. 

Combining similar analyses of all other asymptotic regions we conclude that any trajectory 
of the effective dynamics will always go exponentially to a Kasner trajectory as volume goes 
to zero, with at the most finitely many inter-region traversals (oscillations). 
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VII. SUMMARY AND DISCUSSION 



Bianchi IX models are well known to display the most complex behavior close to a homo- 
geneous singularity but in addition, in the BKL scenario, play a pivotal role in suggesting a 
strongly spatially dependent structure of the spatial slice close to a generic, inhomogeneous 
cosmological singularity. The primary source of this picture is the chaotic behavior of the 
Bianchi IX dynamics which forces infinite fragmentation with no correlations among the 
fragments. Since chaotic behavior itself is a consequence of unbounded growth of the curva- 
ture invariants, i.e. the singularity, and since the loop quantization of diagonal, homogeneous 
models has been shown to be singularity free, it is natural to focus on the modification due 
to the loop quantization of the vacuum Bianchi IX model. 

Since the chaotic dynamics is best understood in the classical setting, we worked with an 
effective classical dynamics from the underlying, singularity free quantum dynamics. This 
was done by exploiting the (always) available continuum approximation together with a 
WKB ansatz for the continuum wave function. The most significant, qualitatively distinct 
effects of the loop quantization got incorporated into the effective potential via the definition 
of inverse triad components together with the quantization ambiguity parameters (j > 
|,0 < / < 1). The effective potential has several qualitative differences from the classical 
potential, (i) Its volume dependence is not in a factorizable form which makes the dynamical 
equations explicitly 'time' dependent, (ii) It does have walls for large volume but 
unlike the classical potential these walls are not confining - they have a finite height and 
a finite extent, (iii) For smaller volume, these walls just disappear, the potential becomes 
negative everywhere and approaches zero in the zero volume limit. Qualitatively one expects 
a simpler, non-oscillatory, behavior asymptotically resembling Kasner behavior since the 
potential vanishes. 

We presented the details of the effective dynamics at two levels. In section El we stayed 
entirely within the domain of validity of the effective model and demonstrated the non- 
asymptotic stability of the Kasner motion to the inclusion of the effective potential. The 
vanishing of inverse triad eigenvalues for small triad eigenvalues was directly responsible 
for this stability. In section IVTl we extrapolated the effective dynamics beyond its natural 
domain of validity into the small volume regime, much as one does in the usual classical 
analysis. We proved in detail that for every initial condition, after at most finitely many 



35 



oscillations, the trajectory will diverge to infinite anisotropics as a Kasner motion. The 
classical case was also contrasted. All the proofs work uniformly for all the allowed values 
of the ambiguity parameters demonstrating the robustness of loop quantum effects. 

To explore the possible implications of the results for the structure of inhomogeneous 
singularities, let us assume the validity of the BKL picture and divide up a spatial slice into 
'cells' which are approximately homogeneous and have (initial) anisotropics parameterized 
by Bianchi IX. We may picture the so decomposed initial spatial shce as a 'sprinkhng of 
points' in the Bianchi IX anisotropy plane, each point corresponding to a cell. The approach 
to an inhomogeneous singularity can now be viewed as evolution of the initial points under 
a Bianchi IX effective (or otherwise) dynamics. 

When the evolution is chaotic (as in the classical case), two nearby points will diverge 
away, effectively loosing initial correlations implied by the smooth geometry of the slice. 
This needs further decomposition and so on. Asymptotically then there would be infinitely 
many points performing un- correlated motion in the anisotropy plane - i.e. a non-interacting 
two dimensional 'gas'. 

When the chaotic motion is replaced by asymptotic Kasner evolution (as with the effective 
dynamics), the points will move in a simpler manner getting distributed into the various 
regions and after some initial crossings between regions they will remain confined to their 
respective regions. Their correlations, just before this stage, will not be washed away and 
we will have an 'interacting gas.' In both cases, all angular sectors are 'filled' suggesting 
that a sufficiently large class of inhomogeneities can be parameterized as a distribution of 
anisotropics to the cells. (The inhomogeneities would be arbitrary in the case of chaotic 
dynamics.) 

Just for contrast assume that Bianchi IX dynamics is governed by the negative of the 
classical potential. Since the walls disappear, this is a trivial way of removing chaos (by 
hand). The stable fixed points of the asymptotic dynamics now include all those for the 
classical potential but in addition the isolated fixed point at a± — 0,a\\ — = — -\/2 is also 
accessible. Now the asymptotic motion compatible with the region would be to the isolated 
fixed point. By permutation, one will have just three possible asymptotic velocities (along 
the dips in the classical potential). All the initially distributed points will cluster in these 
directions and the singularity will look essentially homogeneous. 

We will now use these three scenarios, classical potential, negative classical potential. 
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and effective potential, to shed light on the general structure of classical singularities. In 
the context of the approach to a classical singularity it has been observed j36| that the 
beginning and the end of a universe like ours appear very different from each other. Our 
observable part of the universe starts out very homogeneous and isotropic, but turns much 
more complicated in the future. Black holes form and their singularities make space more 
and more inhomogeneous. Let us now assume that the universe may re-collapse. Then 
by time reversal symmetry the structure of both the initial and the final singularity at 
the boundaries of space-time would have to be the same. If the universe does form black 
holes (which once formed will not disappear in a classical context), then the inhomogeneity 
implied by their presence will not disappear implying a strongly inhomogeneous structure 
for the final singularity. From what we can tell by observations of our own universe, then, 
the situation looks far from time reversal symmetric. Right now space is very homogeneous 
and isotropic at large scales, and is expected to have been even more homogeneous earlier, 
while we do not see a reason for it to homogenize away the structure implied by black holes 
in the future. 

The large degree of homogeneity is usually explained by referring to inflation, which 
expands a small patch of an initial space and thus removes inhomogeneities. (However, 
whether inflation begins at all in an inhomogeneous context is apparently an unresolved 
issue IstI Issj l . Moreover, if the classical BKL dynamics is true, this initial space close to the 
classical singularity would have structure at arbitrarily small scales which would be difficult 
to smooth out even for inflation. If the chaotic behavior leads to a fractal initial state, 
inhomogeneities would not be removed at all due to self-similarity.) If the inhomogeneous 
structure is determined by the chaotic Bianchi IX dynamics, then the primordial anisotropics 
would be completely arbitrary in such a patch. However, there are arguments to the effect 
that any such anisotropics would be washed away before inflation ends js^. This would 
be good for phenomenology in that one does not have to worry about the details prior to 
the beginning of inflation but then one has also lost any contact with whatever quantum 
effects that would be present in the earlier epoch. Alternatively, those effects could be 
influential in determining the beginning of inflation itself. Deciding this may be harder 
when homogeneous patches are un- correlated. In any case, the classical expectation for the 
structure of singularities would be consistent with the expected inhomogeneous end, but 
could be problematic for the beginning. 
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The hypothetical case of a negative classical potential considered above imphes homoge- 
neous structure for the singularity. This would be quite desirable for the initial singularity 
since it is compatible with the observed large scale homogeneity. One may at the most 
have to worry about the isotropization issue. However, it would be in conflict with an in- 
homogeneous final singularity and hence inconsistent with black hole formation during the 
evolution of the universe. Thus merely removing the chaotic Bianchi IX dynamics would 
not be satisfactory. 

The effective dynamics, finally, offers a third possibihty. It admits a sufficiently general 
class of inhomogeneities (thereby accommodating black hole formation) but does not give 
rise to structure on arbitrarily small scales and retains some correlations among the patches. 
This may make starting inflation easier. Alternatively the mechanism of inflation could 
be very different which will not only homogenize but may also retain some imprint of the 
underlying quantum dynamics. The consequences are not worked out at present but it does 
open up other possibilities to be explored. In this aspect, matter will almost certainly play 
a significant role. The influence of matter on the effective dynamics as well as the other 
issues mentioned above will be dealt with elsewhere. 

While the effective dynamics is consistent with our expectations for both the beginning 
and the end of the universe, the apparent time reversal asymmetry remains. This is explained 
by the fact that the situation is, in fact, time asymmetric due to our own position in the 
universe. We can see only some part of it, not the whole space-time, and in particular we 
see only a small part of the beginning. With the current understanding, the observable 
part of our universe can well be part of a classical space-time with a very inhomogeneous 
initial singularity. Since most of the initial singularity is unobservable, however, it is not 
discussed further. The final singularity, on the other hand, is completely unobservable until 
it is reached. If we compare only observable properties from within the universe, we simply 
cannot possibly know enough to tell whether past and future singularities are similar. If we 
compare the theoretical structure of a space-time from outside, then we conclude that in 
fact there is no conceptual difference between the beginning and the end of a generic space- 
time. Only if we compare the observable part of the initial singularity with the theoretical 
expectation for a final singularity does the time asymmetry appear. 
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